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ABSTRACT 



A sample of 198 soft X-ray selected active galactic nuclei (AGN) from the 
■ ROSAT International X-ray Optical Survey (RIXOS), is used to investigate the 

I X-ray luminosity function and its evolution. RIXOS, with a flux limit of 3 x 

i 10~"'^'*erg s~^ cm^^ (0.5 to 2.0 keV), samples a broad range in redshift over 20 deg^ 

0^ ' of sky, and is almost completely identified; it is used in combination with the Ein- 

(— I I stein Extended Medium Sensitivity Survey (EMSS), to give a total sample of over 

600 AGN. We find the evolution of AGN with redshift to be consistent with pure 
I ' luminosity evolution (PLE) models in which the rate of evolution slows markedly 

O I or stops at high redshifts (z > 1.8). We find that this result is not affected by the 

inclusion, or exclusion, of narrow emission line galaxies at low redshift in the RIXOS 
^ ' and EMSS samples, and is insensitive to uncertainties in the conversion between flux 

values measured with ROSAT and Einstein. We confirm, using a model independent 
< Ve/Va > test, that our survey is consistent with no evolution at high redshifts. 

' Key words: surveys - galaxies: active - X-rays: galaxies - cosmology: observations. 



1 INTRODUCTION 

X-ray properties are becoming an increasingly important 
tool for selecting samples of AGN, and sample sizes are 
approaching that of ultraviolet excess (UVX) selected op- 
tical samples. 

The largest currently available sample of serendipi- 
tous X-ray selected AGN comes from the Einstein Ex- 
tended Medium Sensitivity Survey (EMSS, see Stocke 
et al. 1991), with a sky coverage of 778 deg^. How- 
ever, the EMSS is dominated by low redshift objects 
(median z ^ 0.2) due to its high limiting flux (typi- 
cally > 10"^'' erg s"'^ cm~^ in the energy range 0.5 to 
2.0 keV). Deeper Einstein surveys (Primini et al. 1991) 
identified only 11 AGN. With the arrival of ROSAT the 
possibilities for deeper surveys have been realised: the 
Cambridge Cambridge .ROSAT Serendipity Survey covers 
about 4 deg^ and has a flux limit of 2 x 10"^'' erg s~^ cm~'^ 
from 0.5 to 2.0 keV (see Boyle et al. 1995); the survey 
of Boyle et al. (1994) contains 107 broad line AGN and 



probes to flux levels lower than 4 x lO'^^erg s~^ cm~^ 
(0.5 to 2.0 keV), but is only 70% complete at this flux 
limit and covers an area of less than 1.5 deg^ in total. 
This survey is almost devoid of low redshift objects, with 
only 3 AGN of z < 0.4, and median z = 1.5. Even deeper 
surveys, with even smaller sky coverage, (Hasinger et al. 
1993; Branduardi-Raymont et al. 1994) are still in the 
process of optical identification. 

The ROSAT International X-ray Optical Survey 
(hereafter RIXOS, see Mason et al. 1995) occupies a posi- 
tion between the EMSS and the deeper ROSAT surveys, 
with a sky coverage of over 20 deg^ and a limiting flux of 
3 X 10""erg s"^ cm"^ (0.5 to 2.0 keV). It is constructed 
from 81 long (> 8000 seconds) ROSAT pointings made 
with the Position Sensitive Proportional Counter (PSPC) 
at the focus of the X-ray telescope. The RIXOS AGN 
sample consists of 198 objects, the largest ROSAT selected 
AGN sample to date, and is particularly useful for investi- 
gating evolution models because of its broad sampling of 
redshifts. 
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In this paper we present the results of an analysis 
of the differential X-ray luminosity function (XLF) and 
its evolution with redshift, using the RIXOS sample. In 
section 2 we describe briefly the RIXOS survey and present 
details of the RIXOS AGN sample, while in section 3 we 
discuss the log N - log S of RIXOS AGN compared to 
that of the EMSS AGN. In section 4 the methods used in 
this analysis axe explained and in section 5 we present our 
results, which are discussed in section 6. Our conclusions 
are presented in section 7. 

Throughout this paper a Friedmann model universe 
has been assumed, and a value of 50 km s"'^ Mpc~^ has 
been adopted for the Hubble constant Hq; two different 
values for the deceleration parameter, go = and go — 0.5 
have been used. 



2 THE RIXOS SAMPLE OF AGN 
2.1 Observations and Data Reduction 

The X-ray sources in RIXOS were found in a total of 81 
ROSAT PSPC pointings. The target of each observation 
and a small region around it have been excluded from 
the analysis so that RIXOS consists entirely of serendip- 
itously discovered sources. Sources more than 17 arcmin- 
utes offaxis have also been excluded due to their larger 
positional uncertainty and possible masking by the detec- 
tor window support structure. Furthermore, only sources 
detected in the harder ROSAT cncYgj band 0.4 to 2.4 keV 
axe included; the poorer point spread function, interstellar 
absorption, diffuse Galactic X-ray emission, and the in- 
creased contribution of Galactic stars complicate the de- 
tection of extragalactic sources in in the 0.1 to 0.4 keV 
band. Using finding charts from the automatic plate mear- 
suring (APM) facility at the Royal Greenwich Observa- 
tory, Cambridge, optical spectra were taken for all optical 
counterparts within the one- sigrna error circle of each X- 
ray source. If no likely counterpart was found, the optical 
counterparts in the larger two and three-sigma error cir- 
cles were investigated. Where the APM finding charts were 
not sufficient (for example if no optical counterparts ap- 
peared near the X ray source position) CCD images were 
obtained using the Nordic Optical Telescope (NOT), or 
the Isaac Newton Group of telescopes (ING) on La Palma. 
Full details of the optical imaging and spectroscopic ob- 
servations and data reduction are given in Mason et al. 
(1995). 

In this analysis, the term AGN is used to refer to ap- 
proximately the same range of objects as in Maccacaro 
et al. (1991), that is objects with at least one broad ( > 
lOOOkm/s ) emission line and/or [OIII]5007 > [OII]3727. 
Hence the RIXOS AGN sample does include some narrow 
line objects. These criteria have been deliberately chosen 
to avoid any significant difference between the RIXOS and 
EMSS optical selection, allowing the two samples to be 
meaningfully compared and combined. The effect of ex- 
cluding narrow line objects is discussed in section 6. 



Table 1. RIXOS cumulative sky coverage corrected for incom- 
pleteness 

Flux Limit Corrected Optically Number of 
(erg cm~^) Area Identified Fields 
0.5 - 2 keV (deg2) Fraction 



3.0 X IQ-" 


14.16 


93% 


62 


3.5 X IQ-" 


14.36 


95% 


62 


5.0 X 10-" 


14.73 


97% 


62 


6.0 X 10"" 


15.09 


99% 


62 


8.4 X 10-" 


20.04 


99% 


81 



2.2 Construction of the AGN sample used in 
this analysis 

A limiting X-ray flux of 3 x lO^^'^erg s^^ crri^^ (0.5 to 
2.0 keV) was chosen for RIXOS, well above the detection 
threshold for all the J?05ylT fields used. Owing to the con- 
straints of optical telescope scheduling some X-ray sources 
remain unobserved and/or unidentified. In the interests of 
keeping incompleteness and optical selection effects to a 
mimimum, the entire RIXOS survey has not been used. In 
62 of the 81 RIXOS fields, all objects have been identified 
or observed speetroseopieally to the intended flux limit of 
3.0 X lO^^^crg s"^ cm-2 (0.5 to 2.0 keV), wfiile in each of 
the other 19 fields some, but not all, of the X-ray sources 
have been observed to this limit; these 19 fields are, how- 
ever, fully observed to a flux of 8.4 x IQ-^'^erg s~^ cm"^ 
(0.5 to 2.0 keV) and are included in the RIXOS AGN 
sample with this flux limit. Overall spectroscopic com- 
pleteness of the fields used in the RIXOS AGN sample 
is 93%; the remaining 7% of sources which are unidenti- 
fied are those for which the optical counterpart (s) were 
too faint for us to obtain reliable optical spectra. We 
have made the assumption that the fraction of uniden- 
tified sources which axe AGN is the same as that for the 
identified sources. Accordingly, the sky area used for this 
analysis has been corrected, in a similar fashion to that of 
Boyle et al. (1994), by multiplying the area by the fraction 
of sources identified; as the unidentified fraction is small, 
this has only a small effect on our results. Since optical 
completeness is a function of flux limit we have calcu- 
lated the effective sky area at 3.0 x 10"^'' erg s~^ cm"^, 
8.4 X 10~^'*erg s"^ cm"^, and three intermediate fluxes 
corresponding to significant changes in spectroscopic com- 
pleteness. Again, the high level of completeness in RIXOS 
makes this a small correction, which has only a small ef- 
fect on our results. The number of i? 05^4 T fields, corrected 
sky coverages and identified fractions at their respective 
limiting fluxes axe listed in Table 1. 

The redshift distribution, N{z), of the RIXOS AGN 
sample is shown in Fig. 1. The sample has a significantly 
higher median redshift, 0.6, than the EMSS (0.2). 

To obtain the largest possible working sample of 
AGN, the RIXOS and EMSS surveys have been combined 
coherently (Avni & Bancall 1980) to give a total of over 
600 AGN. This combined sample will be referred to as 
'RIXOS + EMSS' hereafter. To correct the EMSS sam- 
ple for incompleteness, the EMSS 'expected' AGN (see 
Maccacaro et al. 1991) have been included in the EMSS 
and EMSS -|- RIXOS samples. Throughout this analysis a 
power law X-ray spectrum has been assumed, fv oc 1/'"^^ . 
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Figure 2. X-ray Luminosity and redshift of RIXOS AGN 
(closed squares) and EMSS AGN (open triangles). Luminosities 
have been calculated using go = 0. 



Figure 3. Integral log N - log S of RIXOS AGN (closed 
squares) and EMSS AGN (open triangles). 

For comparison of the results presented in this paper with 
those of Boyle et al. (1994) and Maccacaxo et al. (1991), 
the analysis has been performed using ax = 1- The first 
results from X-ray colour analysis (see Mittaz et al. 1995) 
indicate that this is quite representative of the RIXOS 
AGN, which have a median ax ~ 1-01. The scatter of ax 
around the value 1.0 appears to be larger in the RIXOS 
sample than was found by Macaccaro et al. (1988) for 
EMSS AGN, although the hardening of source spectra 
with redshift (sec Francis 1993) does not appear to be sig- 
nificant in the RIXOS sample and has not been included 
in this analysis. The X-ray luminosity - redshift {Lx,z) 
distribution for the RIXOS and EMSS AGN are compared 
in Fig. 2. As expected from a deeper survey, the RIXOS 
AGN typically have lower luminosity and/or higher red- 
shift than the EMSS AGN. 



3 LOG N- LOG S 

The survey of Boyle et al. (1994) contains a larger num- 
ber of AGN than would be predicted by extrapolating the 
EMSS log N - log S to lower fluxes; these authors sug- 
gested that this may be partly due to an error in their con- 
version from ROSAT to Einstein fluxes. In practice there 
are a number of factors which might affect the relative 
source counts and/or fiuxes of the two samples, including 
differences in the source detection and parameterisation 
algorithms, uncertainties in the detector response matri- 
ces used to transform count rates to fluxes, uncertainties 
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in the spectral form used to convert from fluxes in the Em- 
stem 0.3 to 3.5 keV band to the ROSAT 0.5 to 2.0 keV 
band and the effects of incompleteness or other biases. 

To assess these effects Figure 3 compares the integral 
log N-\og S for the RIXOS and EMSS AGN. This plot 
is derived using a conversion factor (CF) of 1.8 between 
the fluxes in the two bands which is appropriate for a 
powerlaw spectrum of aa, = 1 if we use the standard pub- 
lished response matrices for ROSAT and Einstein. It is 
clear that RIXOS has a larger number of sources than the 
EMSS by about 30% at almost all fluxes when using this 
conversion. Although RIXOS samples a larger area than 
the EMSS at low fluxes, there is significant overlap in the 
fluxes of objects found in the two samples, and there is 
thus no physical reason why the log N - log S relations of 
the two surveys should be different. 

Further investigation of the reeisons for this discrep- 
ancy is beyond the scope of this paper. Instead wo adopt 
an approach whereby we parameterise the discrepancy em- 
pirically, and investigate to what extent the uncertainty in 
this number affects the results on AGN evolution when we 
combine the RIXOS and EMSS samples. The device we 
use for this empirical parametcrisation is to find the CF 
for which the two log A'^- log S relations arc consistent. At 
the fiux fimits of the EMSS and RIXOS the N(S) relation 
is well fit by a power law. Using maximum likelihood, the 
two samples have been fitted simultaneously with a single 
power law slope but different normalisations. 
rIN 

— = ksS"' EMSS objects, S(0.3 - 3.5keV) 



UrS"^ RIXOS objects, S(0.5 - 2.0keV) 



dN 
dS 

There are two free parameters in this fit, the power law 
slope 7 and the difference between the two normalisations 
(^B/fefl); the actual normalisations are found by requiring 
that the number of objects predicted by the log N log 
S for the sky coverage of RIXOS plus the sky coverage of 
the EMSS is equal to the total in RIXOS plus the total in 
the EMSS. The normalisation difference should be related 
simply to the empirical CF from ROSAT to Einstein fluxes 
by 

CF = (kE/kn)'^^''-'^ 

We assumed upper flux limits of 2 x lO^^^crg s^^ cm~^ 
(0.3 to 3.5 keV) and 10""crg s"^ cm"^ (0.5 to 2.0 keV) 
for the EMSS and RIXOS respectively, to reflect the selec- 
tion against very bright sources in these surveys. This has 
only a small effect on the results; changing the upper fiux 
limits to any reasonable value (or removing them) changes 
the best fit by only a fraction of the 1 a statistical errors 
quoted below. The best fit slope 7 is found to be 2.56, con- 
sistent with the slopes found when the two samples are fit 
independently, 2.61±0.06 for the EMSS and 2.45±0.11 for 
RIXOS, where errors are 1 a. Since the slope of the EMSS 
log N- log S curve is actually steeper than (although con- 
sistent with) that of RIXOS, there is no evidence from this 
comparison to support previous claims about incomplete- 
ness in the EMSS (e.g. Franceschini et al. 1994) in which 
incompleteness is thought to be a problem at low fiuxes. 
The best fit CF is found to be 1.47 ±0.11; again errors axe 
1 a. The standard conversion factor of 1.8 is thus signifi- 
cantly different (rejected at >99.0%) from that found by 
matching the log N - log S curves. There is, however, no 



evidence from the RIXOS and EMSS log N - log S that 
the CF should be as small as 1.0, a possibility considered 
in Boyle et al. (1994). 

To assess the impact of the different conversion factors 
on AGN evolutionary models, values of 1.47 and 1.8 have 
both been used in the subsequent analysis and the results 
compared. Note that all fluxes and luminosities quoted in 
this paper are 0.5 to 2.0 keV, i.e. the Einstein fluxes have 
been converted to the ROSAT flux band. 



4 THE LUMINOSITY FUNCTION 

The techniques that we have applied to the data are dis- 
cussed in this section, while the results of the analysis are 
described in section 5. 

A flux limited survey samples different luminosity 
ranges at different redshifts. The AGN population as a 
function of luminosity and redshift is represented by the 
XLF, 4>{L,z), which is defined as the number of objects 
detected per unit comoving volume per unit luminosity 
interval, i.e. 



<t){L,z) 



-iL,z) 



dVdL' 

Evolution in the space density or luminosity of AGN 
with redshift is seen as a change in the XLF. 

The two simplest forms for evolution are pure den- 
sity and pure luminosity evolution (PLE) . In pure density 
evolution, the number of AGN per unit comoving volume 
is assumed to evolve while the distribution of luminosities 
remains constant. In PLE, the number of AGN per unit 
comoving volume remains constant, while the luminosity 
of each AGN evolves. 

Pure density evolution models do not predict the flat- 
tening at low fluxes which is seen in the best log N - log 
S curves currently available, in both the optical (Boyle et 
al. 1988) and X-ray bands (Branduardi - Raymont et al. 
1994, Hasinger et al. 1993). For this reason, recent evolu- 
tionary models have been based on luminosity evolution of 
some form. Only PLE models are considered in this paper. 

As shown by Maccacaro et al. (1991), the XLF can 
be modeled as two power laws with a break luminosity, 

Lbreak , SUCh that 

<j) — KlL~^^ L < Lbreak 

4> = K2L~^^ L > Lbreak 

where Ki and K2 are normalisations of the two power 
laws. Since we require the luminosity function to be con- 
tinuous, the two normalisations are not independent. A 
single normalisation Ki is adequate, since 

In the framework of a PLE model the XLF retains 
its shape at all redshifts, hence the XLF at any redshift 

depends on only the evolution law and the XLF at zero 
redshift (hereafter z = XLF) . 



4.1 Maximum likelihood and 2D K.S. testing 

The maximum likelihood method (Marshall 1984) has 
been used to obtain best fit evolution parameter (s) and 
z = XLF for each evolution model, utilising the full 
RIXOS + EMSS sample of AGN. This technique involves 
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simultaneously fitting the evolution and the z = G XLF. 
There are four or five free parameters (three from the z 
= XLF, 71, 72, and Lbreak, plus either one or two from 
the evolution model) in fitting the models considered in 
this paper. The normalisation of the z — XLF is set so 
that the total number of objects predicted by the model 
is equal to the number in the sample, and is not a free 
parameter in the fit. 

To test the acceptability of evolutionary models, the 
two dimensional Kolmogorov-Smirnof (2D K.S.) test has 
been used. Two alternatives have been used in the lit- 
erature, the test described by Peacock (1983), and the 
test of Fasano & Franceschini (1987), which is also de- 
scribed in Press et al. (1992). It is important to use the 
test most effective at rejecting poor models and accepting 
good models for the RIXOS and EMSS datasets, to ensure 
confidence in the results. To assess the most appropriate 
test to use, both were applied to simulated data as follows: 

The two best fit XLFs and evolution models (see 
section 5.1) from Maccacaro et al. (1991) were used to 
produce 200 simulated samples of AGN, 100 from each 
evolution model. Each sample was constructed using two 
completeness limits, appropriate to the RIXOS and EMSS 
samples respectively. The number of objects in each sim- 
ulated sample was between 400 and 600, similar to the 
number of objects in the combined RIXOS + EMSS sam- 
ple, and the tests were performed over exactly the same 
plane as those used on the real sample (see below) . In this 
way, the 2D K.S. tests were evaluated under very simi- 
lar conditions to those under which they were actually to 
be applied. Each sample was tested once against its own 
parent XLF and evolutionary model, and once against the 
other, so that 200 tests were performed with correct mod- 
els and 200 with incorrect models. The Peacock test re- 
jected 100 incorrect models (50%) at the 95% level and 6 
correct models, while the Fasano & Franceschini test re- 
jected 84 incorrect models (42%) and 10 correct models. 
This indicates that the 2D K.S. test of Peacock is more 
efficient at distinguising between good and bad models of 
our data. Because the EMSS and RIXOS contain a very 
large range of sky coverage at different completeness lim- 
its, the combined RIXOS -|- EMSS sample has a slightly 
lower redshift - de-evolved luminosity correlation coeffi- 
cient (typically 0.43, but dependent on cosmological and 
evolutionary model) than the simulation models (typically 
0.5 to 0.7). According to Fasano & Franceschini (1987), 
their test reaches maximum efficiency when compared to 
that of Peacock (1983) at higher correlation coefficients, 
and so it is reasonable to assume the results obtained 
from our simulations should hold for the actual data to 
be tested. The test of Peacock (1983) has therefore been 
used in this paper. 

Models have been tested in the redshift and de- 
evolved luminosity plane over a range that includes all 
parameter values found in RIXOS and the EMSS, (0 < 
2 < 3.5, 10*° < Lo < W*'^) where Lo is the de-evolved 0.5 
- 2.0 kcV luminosity in erg s^^; note that Lo is model de- 
pendent. Since there is no selection criterion in this anal- 
ysis, based on observed X-ray luminosity, which would 
correspond to the Mb < — 23 requirement often used in 
optical QSO surveys, the test has not been performed over 
an interval in observed luminosity. Imposing a lower limit 



to observed luminosity in this way would introduce im- 
plicit density evolution to a PLE model (see Kassiola and 
Mathez 1990), and is hence undesirable. 

4.2 1/Va and < Ve/Va > 

An estimate of the behaviour of the XLF can be gained in 
a model-independent way using the 1/Va statistic (Avni 
& Bahcall 1980) and plotting the XLF in distinct red- 
shift bins, as in Maccacaro et al. (1991). Of course, any 
evolution occuring within these distinct bins will not be 
apparent in this type of treatment. 

A more quantitative treatment can be made using 
the < Ve/Va > statistic (Avni & Bahcall 1980). To in- 
vestigate single parameter PLE models, without simulta- 
neously modelling the z = XLF, we have carried out 
< Ve/Va > tests in distinct redshift shells. In this case the 
test is used to obtain acceptable evolution parameters for 
the given model, and the bins have been chosen so that 
within each shell the evolution parameter has a 68% con- 
fidence region of about ±20% its value over the entire red- 
shift range; this is a good compromise between resolution 
in redshift and constraint of the evolutionary properties 
of each bin. The practical aspects of the < Ve/Va > test 
within redshift shells are described in Delia Ceca et al. 
(1992). 

We have also used the < Ve/Va > test to examine 
evolution at high redshift in a model independent way. 

In this case the < Ve/Va > test has boon applied in the 
redshift interval Zb < z < 3.5 where zt is varied. Here, the 
test is used with no evolutionary model, and is capable 
of determining whether the data are consistent with the 
no evolution hypothesis between Zb < z < 3.5, and if 
not, whether the luminosity function is increasing with 
redshift (< Ve/Va > > 0.5), or decreasing with redshift 
(< Ve/Va > < 0.5). This use of the < Ve/Va > test is 
described in Dunlop and Peacock (1990). 



5 RESULTS 

The best fit evolution parameters and luminosity func- 
tions have been obtained from the combined RIXOS + 
EMSS sample. As a check for consistency between the two 
samples they have been tested for goodness of fit both in- 
dividually and in combination. Table 2 shows the results 
for the maximum likelihood and 2D K.S. tests applied 
to the models that have been investigated. Errors quoted 
were obtained using the method of Lampton, Margon and 
Bowyer (1976), and correspond to Ax^ = 1, i.e. 68% con- 
fidence intervals for one interesting parameter. 

5.1 The simplest models 

The two simple PLE models used by Maccacaro et al. 
(1991), are power law evolution 

L = LoX (l + zf 

and evolution which is exponential with look back time 

L = LoX e^^ 

where C is the evolution parameter and t is look back 
time. 
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Table 2. Results of fitting evolution models 

model CF go z,^t C Ci Xi| 71 72 logw Ixrb^ P(>D) P(>D) P(>D) 

(^-breafe)! RIXOS RIXOS EMSS 
+EMSS 



(l + zf 


1.8 


0.0 


— 


2.66 


— 


2.18 


1, 


.68 


3.38 


3.59 


8.24 


0.020 


0.26 


0.036 




1.8 


0.0 




4.74 




1.59 


1, 


.64 


3.23 


3.29 


7.31 


0.033 


0.047 


0.033 




1.8 


0.0 


1.82 


2.91 




1.83 


1, 


.65 


3.30 


3.49 


6.88 


0.27 


0.62 


0.096 




-L.O 


n n 
u.u 




1 1 n 




i.OU 




.uz 


97 


'\ ^9 


u.^ / 


u.ou 




n 1 7 


(1 + 2)C 


1.8 


0.5 


- 


2.35 


- 


1.75 


1, 


.62 


3.38 


3.57 


4.56 


0.015 


0.62 


0.017 




1.8 


0.5 




3.82 




1.37 


1, 


.63 


3.32 


3.20 


4.71 


0.019 


0.21 


0.0048 


{l + zf 


1.8 


0.5 


1.42 


2.94 


— 


1.22 


1, 


.57 


3.30 


3.38 


4.18 


0.69 


0.24 


0.25 




1.8 


0.5 


- 


1.07 


-0.238 


1.48 


1, 


,60 


3.33 


3.46 


4.20 


0.58 


0.58 


0.21 


(l + zf 


1.47 


0.0 




2.46 




2.00 


1, 


.64 


3.23 


3.64 


8.28 


0.055 


0.55 


0.042 




1.47 


0.0 




4.44 




1.81 


1, 


.65 


3.16 


3.40 


7.71 


0.016 


0.095 


0.064 


(1 + z)^ 


1.47 


0.0 


1.82 


2.72 




2.19 


1, 


.66 


3.23 


3.60 


7.35 


0.22 


0.86 


0.19 


IqCz + Ciz'^ 


1.47 


0.0 




1.04 


-0.219 


1.92 


1, 


.64 


3.20 


3.58 


7.03 


0.24 


0.65 


0.20 


il + z)C 


1.47 


0.5 




2.19 




1.64 


1, 


.59 


3.23 


3.62 


4.88 


0.049 


1.2 


0.019 




1.47 


0.5 




3.47 




1.52 


1, 


.61 


3.19 


3.31 


4.86 


0.026 


0.24 


0.045 


{\ + z)C 


1.47 


0.5 


1.41 


2.73 




1.39 


1, 


.57 


3.19 


3.48 


4.57 


0.22 


0.38 


0.14 




1.47 


0.5 




1.04 


-0.238 


1.54 


1, 


.58 


3.22 


3.51 


4.69 


0.38 


0.80 


0.18 


itERRORS 






























(1 + ^)'^ 








0.08 






0.08 
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T = l-l/{l + zf^^ 90=0.5 

The 2D K.S. test rejects both these simple models at the 
95% level for the combined sample except the go = 0, CF 
= 1.47 case where the power law model is just acceptable 
at the 95% level. Note that the case where the 2D K.S. 
probability P(> D) is 1.2 indicates that the difference 
between the expected and actual distributions is small, 
not that it is zero; it is possible for the test of Peacock 
(1983) to produce values for the P(> D) larger than 1, in 
which case it cannot be regarded as a probability, although 
the implication that a model is a good fit if P(> D) is high 
is certainly true. 

Figs 4 and 5 show binned 1/Va luminosity functions 
in redshift shells against the best fit power law and expo- 
nential evolution model luminosity functions. Figs 6 and 
7 shows evolution parameters acceptable at 68% for these 
two evolution models derived from the < Ve/Va > test in 
redshift intervals. These figures suggest that the exponen- 
tial model fails because it requires unacceptably rapid evo- 
lution at low redshift, (in Fig. 7 the low redshift evolution 
parameter lies below any value which would be consistent 
with z > 0.4), while the power law model fails because it 
overpredicts evolution at high redshift, (for z > 1.5 the 
model curves in Fig. 4 lie well above the data). Evolution 
is slower at high redshift in the exponential model than in 



the power law model, and hence in Fig. 5 the exponential 
model appears less discrepant ai z > 1.5 than the power 
law model in Fig. 4. Figs 4 to 7 have been constructed 
using qo = 0, CF = 1.8. 



5.2 Two Parameter Evolution Models 

As both single parameter models are rejected by the 2D 
K.S. test, more complex models have been investigated. 
From Figs 4 and 5 it appears there is little difference in 
the slope of the luminosity function in different redshift 
bins, so models in which the luminosity function changes 
slope with redshift have not been considered. A power law 
model with a redshift cutoff, Zcut, where 
L = Lo X (1 + z)^ z < Zcut 

L = Lo X (1 + Zcut)^' Z > Zcut 

and alternatively, a polynomial evolution of the form 
L = Lo X iQ^c^+c^^') 

both have two free parameters for the evolution (plus three 
for the luminosity function). Both of these models are ac- 
cepted at the 95% level by the 2D K.S. test for the com- 
bined RIXOS -I- EMSS samples for both values of qo and 
CF; there is little justification, from the 2D K.S. probabil- 
ities obtained, to prefer one of the models. The similarity 
of the two models for z < 1.5 is illustrated by the pre- 
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Figure 6. Power law evolution parameter C in redshift bins 
using the < Ve/Va > test for qo = and CF = 1.8. 



Figure 4. Binned < 1/Va > XLF of the RIXOS + EMSS 
sample and best fit power law evolution model XLF (dashed 
lines) for qo = and CF = 1.8. 
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Figure 5. Binned < 1/Va > XLF of the RIXOS + EMSS 
sample and best fit exponential evolution model XLF (dashed 
lines) for go = and CF = 1.8. 
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Figure 7. Exponential evolution parameter C in redshift bins 
using the < Ve/Va > test for qo = and CF = 1.8. 



dieted N{z) relations plotted in Fig. 1. The two evolution 
models are however radically different in shape beyond 
this redshift; in the polynomial evolution model luminos- 
ity declines after 2, and at z = 3.5 the expected number 
of objects differs by a factor of 5 for the two models. A 
larger sample of z > 2 objects would have the potential to 
discriminate between the two models. 



5.3 Evolution at high redshift 

It is seen in the previous section that PLE models in which 
evolution ceases or changes direction at high redshift are 
found more acceptable than models where evolution con- 
tinues, for both values of go and CF, indicating that evo- 
lution at « > 2 must be absent or slow compared to that 
at low redshift. As further evidence, the results of the 
< Ve/Va > test from z = zt to z = 3.5 for the combined 
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Figure 8. < Ve/Va > test in the redshift interval z ■ 
= 3.5 using the RIXOS + EMSS sample for qo = 0. 



Zb to 2 Figure 9. < Ve/Va > test in the redshift interval z = 
= 3.5 using the RIXOS + EMSS sample for go = 0.5. 



Zb to z 



RIXOS + EMSS sample arc shown in Figs 8 and 9. Both 
were constructed using a CF of 1.8; when 1.47 is used, 
< Ve/Va > is a few percent lower; note that above Zb = 
2.4 there are only 4 objects included in the test. From Zb = 
1.67, with 30 AGN, the < K/K > test shows the data to 
be consistent at the 68% level with no evolution for both 
values of qo and CF. The < Ve/Va > test used in this way 
is model independent, and this result is not restricted to 
PLE models. 



6 DISCUSSION 

Although changing the ROSAT to Emstem CF from 1.8 to 
1.47 does not appear to affect the choice of model, (i.e. the 
power law with evolution cutoff and polynomial models are 
both significantly more acceptable than the simple power 
law and exponential models), it does have a significant 
effect on the best fit evolution paramcter(s) and the z = 
XLF. In all cases the best fit parameters for CF = 1.47 
are outside the 90% confidence region of the parameters 
for CF = 1.8; evolution is slower, and the steep part of the 
XLF is less steep if the CF of 1.47 is used. The slope of 
the low luminosity region of the XLF is affected little by 
the choice of CF. 

The evolution parameter, C, for the exponential 
model is consistent at 68% with the EMSS value (C = 
4.18±0.35 for go = 0) only for CF = 1.47, and its low K.S. 
probability confirms the result of Delia Ceca et al. (1991) 
that this model is a poor description of AGN evolution. 
For the power law model, the value of C obtained here 
is consistent with that found by Maccacaro et al. (1990) 
(2.56 ± 0.17 for qo = 0) for both values of CF. The in- 
troduction of a cutoff in evolution aX, z = 1.8 significantly 
improves the 2D K.S. probability of the power law model 
for the RIXOS + EMSS sample; the last column of Table 
2 shows that the EMSS data alone is also better fit with 



the evolution cutoff. Including an evolution cutoff at z = 
1.8 in a < Ve/Va > test to the EMSS data with a power 
law evolution model and qo = gives C = 2.74 ± 0.20, al- 
most identical to the RIXOS+EMSS value for CF = 1.47 
and still consistent for CF — 1.8. 

Comparing the CF — 1.8, go =0 evolution parame- 
ters and « = XLFs from the RIXOS+EMSS sample with 
those of Boyle et al. (1994), we find similar values for all 
but the exponential evolution model which evolves faster 
in Boyle et al. (1994). For qo = 0.5, there is less agreement, 
with the 2 = XLFs of the exponential, power law, and 
polynomial evolution laws having significantly lower val- 
ues for 7i (i.e. flatter slopes at low luminosity) in Boyle et 
al. (1994). The best fit go = 0.5 power law with evolution 
cutoff model in Boyle et al. (1994), whilst having a similar 
z = XLF and evolution rate C, has a much higher cutoff 
redshift {zcut = 1.7), than the value found here, Zcut = 1.4 
for go — 0.5. 

It is notable that almost all of the models considered 
in this paper are found more acceptable to the 2D K.S. 
test than in Boyle et al. (1994), in which almost all mod- 
els are rejected at > 99%. Boyle et al. used the test of 
Fasano & Franceschini (1987), while we have used that 
of Peacock (1983). However we have also tested our mod- 
els against the RIXOS -|- EMSS dataset using the test of 
Fasano & Franceschini, and do not find them to be re- 
jected at > 99%. In Boyle et al. (1994), the inclusion of 
narrow line objects gave significantly increased 2D K.S. 
probabilities, (i.e. a better fit). However, narrow line ob- 
jects which would be classified by Stocke et al. as AGN 
have been included throughout this analysis. If we exclude 
these objects, which occurr in both the EMSS and RIXOS, 
we do find lower 2D K.S. probabilities, although the poly- 
nomial and power law with evolution cutoff models are 
still acceptable at the 95% level while the simple power 
law and exponential evolution models are not. For the z = 
XLF, typically 71 is reduced by about 0.1 and Lbreak is 
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increased by 25% to 50% depending on the specific PLE 
model and choice of cosmology; the best fit values of 72 
and Zcut change by no more than 0.02 and the evolution 
parameters remain within the errors quoted in table 2. 
There arc no broad lino objects with Lo < 2 x 10*^ erg s~^. 

The poor model fits of Boyle et al. (1994) may be 
attributable to the conversion between ROSAT and Ein- 
stein measured fluxes; the ROSAT sample of Boyle et al. 
(1994) combined with the EMSS AGN using a OF of 1.8 
shows an evolution rate too high to be consistent with the 
EMSS sample alone. As we have discussed in section 3, 
the effective OF from ROSAT to EINSTEIN fluxes may 
be significantly lower than 1.8, and as seen in table 2 a 
lower CP has the potential to reduce the evolution rate of 
combined ROSAT and EINSTEIN data, and hence could 
improve the self consistency and model fits of Boyle et al. 
(1994). 

The log A'^ - log S at faint fluxes derived from the 
CP =1.8 models tested above, and extended to 2 = 4 
are shown in Pig. 10. None of the model curves exceed 
the total faint X-ray log N - log S obtained by fluctua- 
tion analysis in Hasinger et al. (1993), or Barcons et al. 
(1994). Notably, the log A'^- log 5 curves at faint fluxes for 
PLE are separated strongly by the value of go used, while 
the specific choice of PLE model has a comparatively mi- 
nor effect. In a go = 0.5 universe, AGN undergoing PLE 
should represent between 30% and 45% of all sources with 
S > lO^^'^crg s~^ crn^^; in contrast, in a go = uni- 
verse AGN evolving in this way would be expected to 
constitute between 55% and 100% of these sources. Above 
3 X 10~^''erg/s/cm^ the log N-log S curves from the PLE 
models are all very similar and represent the data well. 

The AGN contribution to the 1-2 keV X-ray back- 
ground has been calculated for < « < 4, 10'"' < Lo < 
10^^ (where Lo is the de-evolved 0.5 - 2 keV luminosity in 
erg s~^) for all the models tested in section 5 and is shown 
in table 2 (column entitled Ixrb)- As expected from the 
log TV - log S predictions, the contribution of AGN to the 
X-ray background from these models has a stronger de- 
pendence on the value of go than the choice of PLE model. 
The values for the AGN X-ray background intensity given 
in table 2 are in good agreement with those of Boyle et al. 
(1994). A recent measurement of the X-ray background 
(Chen, Fabian and Gendreau 1996) using both ASCA and 
ROSAT iowa^d an intensity of 1.46 xlO~*ergs"^cm"^sr"^ 
(1-2 keV). Using this value, our acceptable models (power 
law with evolution cutoff and polynomial) predict that 
AGN account for between 44% and 50% (go = 0) or 29% 
and 32% (go = 0.5) of the 1-2 keV X-ray background, 
where these ranges include the uncertainty in CP. Taking 
the value of 1.25 x 10~*ergs"^cm"^sr"^ for the 1-2 keV 
X-ray background (Hasinger 1992) used by Boyle et al. 
(1994), the contribution from AGN rises to between 52% 
and 60% (go = 0) or 33% and 38% (go = 0.5). 



7 CONCLUSIONS 

We have investigated the evolution of the XLF with 
redshift using a new sample of 198 X-ray selected 
AGN with a spectroscopic completeness of 93% at 3 x 
10-"erg s"^ cm-2 (0.5 to 2.0 keV). We find PLE mod- 
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Figure 10. Log N - log S of the evolution models extrapo- 
lated to faint fluxes. The thick lines represent the la confl- 

dence limits to the total log A'' - log S of all sources (i.e. not 
only AGN) at faint fluxes determined by fluctuation analysis 
of a deep ROSAT ?Le\d by Barcons et al. (1994). 

els, consistent with our data, in which the XLF declines 
or ceases to evolve beyond z ~ 1.8, and wc find no evi- 
dence for evolution beyond this redshift from the model 
independent < Ve/Va > test. We have shown that nar- 
row emission line galaxies at low redshift contaminating 
the AGN sample have little effect on the derived evolution 
properties. Furthermore, these conclusions are insensitive 
to the uncertainty in conversion of fiux from the Einstein 
to iiOyS^Tpassbands, although this has a significant effect 
on the z = Q XLF and evolution rate for z < 1.8. 
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